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Abstract 

The classical solutions to higher dimensional Yang-Mills (YM) systems, which form part of 
higher dimensional Einstein- YM (EYM) systems, are studied. These are the gravity decoupling 
limits of the fully gravitating EYM solutions. In odd spacetime dimensions, depending on the 
choice of gauge group and its representation, these are either topologically stable or are unstable. 
Both cases are analysed, the latter numerically only. In even spacetime dimensions they are 
always unstable, describing saddle points of the energy, and can be described as sphalerons. This 
instability is analysed by constructing the noncontractible loops and calculating the Chern- 
Simons (CS) charges, and also perturbatively by numerically constructing the negative modes. 
This study is restricted to the simplest YM system in spacetime dimensions d = 6, 7, 8, which 
captures the qualitative features of the generic case. 



1 Introduction 



Gravitational and non Abelian gauge fields occur in low energy effective actions [1] of superstring 
theory and supergravities. On the other hand classical solutions to the Einstein- Yang-Mills (EYM) 
system, especially black holes, have an important role to play in quantum gravity [2]. These effective 
actions consist, in addition to the usual Einstein-Hilbert and Yang-Mills (YM) systems, of higher 
order terms in both the gravitational curvature, and the YM curvature and its gauge covariant 
derivatives. Moreover some of the theories in which such terms are present are defined in higher 
dimensional spacetimes, namely on D-branes [3]. Thus, the study of classical solutions to higher 
curvature EYM models in higher dimensions is of physical interest. 

It turns out that in higher (than 3 + 1) dimensions, the higher order YM curvature terms play a 
much more important role than do the corresponding gravitational terms, e.g. Gauss-Bonnet terms. 
Inclusion of the latter does not seem to alter the qualitative properties of the classical solutions, while 
the absence of higher order YM terms prevents the existence of such solutions due to the Derrick 
scaling requirements not being satisfied. It is for this reason that we restrict our considerations in 
this paper to the inclusion of higher order YM 1 only. Our aim being the study of the stability (and 
its absence) in these systems, we ignore the gravitational terms entirely since the most important 
mechanism affecting stability/instability is characterised by the YM sector alone. 

Concerning the higher order YM curvature terms in the string theory effective action, the situation 
is complex and as yet unresolved. While YM terms up to F 4 arise from (the non Abelian version 
of) the Born-Infeld action [5], it appears that this approach does not yield all the F 6 terms [6]. 
Terms of order F 6 and higher can also be obtained by employing the constraints of (maximal) 
supersymmetry [7]. The results of the various approaches are not identical. 

Given the evolving stage which higher order curvature YM terms are in, and motivated by the 
technical requirements for the construction of classical solutions, we have restricted our considerations 
to one particular family of higher curvature YM systems. The criterion is that only the second 
power (and no higher power) of the velocity field doAi occur in the Lagrangian. This constrains our 
choice to one where the coefficients of each F 2n term is fixed by the requirement that only totally 
antisymmetrised curvature n-forms are employed. We have no physical justification for this, but we 
hope that the ensuing qualitative results hold also in the more general, and as yet not definitely 
fixed, cases. In this family of YM systems the number of higher order terms that can arise in any 
given dimension is limited due to the imposed antisymmetry. 

The other constraint we apply is that in any given dimension, we truncate the series of F 2n terms, 
to the minimum required to satisfy the Derrick scaling requirement. Our justification here is that 
higher order terms become more important at high energies so in the low energy effective action it 
is sufficient to keep the lowest order terms. 

It is well known that there exist static regular [8] and black hole [9, 10] solutions to the Einstein- 
Yang-Mills (EYM) system in d = 4 spacetime dimensions. To construct static solutions to gravitating 
Yang-Mills (YM) systems in spacetime dimensions d > 5 [11, 12], terms of higher order in the YM 
curvature must be included. Such terms appear in various low energy effective actions of string 
theory. For our practical purposes, a particulary useful family of such YM systems is 

^g^™ 2 (i) 

1 Higher order gauge field curvatures arise also as quantum corrections, see e.g. [4], but this is not the source of 
higher order terms we have in mind here. 
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which is the sum of the terms in the YM hierarchy [13], with F(2p) being the p fold totally antisym- 
metrised products of the YM curvature, F(2) = F^ u , in this notation. Clearly, the highest value P 
of p in (1) is finite and depends on the dimensionality d of the spacetime. In [11, 12] we chose the 
simplest possibility P = 2, since we restricted our study to 8 > d > 5. The study of the classical 
solutions to the systems (1), which turn out to be the gravity-decoupling solutions of EYM systems 
in these dimensions, is the aim of this paper. This is an important aspect of the study of higher 
dimensional EYM solutions. 

To complete the definition of the models (1) the gauge group G must be specified. With the aim 
of constructing static spherically symmetric solutions in d spacetime dimensions, the smallest such 
gauge group is G = SO(d — 1) 2 . Static finite energy solutions to the systems (1) may or may not 
be topologically stable depending on the dimensionality d of the spacetime and the choice of gauge 
group G. 

The purpose of the present work is to resolve the question of stability or instability of these 
spherically symmetric solutions quantitatively. These are the gravity decoupling limits of the static 
and regular solutions to EYM systems in spacetime dimensions d > 6 [11]. In the d — 5 case [12], 
as in the usual d = 4 case [8], there are no gravity decoupling solutions, unlike in the case of the 
EYM-Higgs (EYMH) systems in d — 4 whose regular solutions [14, 15] tend to the 't Hooft-Polyakov 
monopole in the flat space limit. In both cases, higher dimensional EYM and EYMH in d = 4, the 
models feature an additional dimensional constant. As a result gravitating solutions exist only for a 
finite range of values of the gravitational constant. Here, we restrict to the non-gravitating models 
only. 

For simplicity, we restrict to the model considered in [11], (1) with p = 1 and 2 only, in spacetime 
dimensions d = 6, 7 and 8. Solutions to this model on d = 6 Euclidean space were constructed in [16] 
in a different context. Consistently with our requirement for the imposition of spherical symmetry 
we choose G = SO(d — 1) and G = SO(d), and because of the higher order of nonlinearity in the 
p > 1 models the choice of the concrete representations of G is important. 

Except for the case where there exists a nonvanishing Chern-Pontryagin (CP) charge in the 
spacelike dimensions, these solutions are unstable. In the particular cases where G and its represen- 
tation is chosen such that the Chern-Simons (CS) charge is nonzero, this instability will be that of 
a sphaleron [17, 18]. The quantitative study [19] of the latter will be the major task below. 

In section 2, we present the detailed models and the corresponding CS densities. In section 3, 
spherical symmetry in the d — 1 spacelike dimensions is imposed and the residual one dimensional 
models are displayed in 3.1 where the gauge invariance of the Ansatz is also stated. The static 
equations are given in 3.2 and the corresponding CS densities in 3.3. In 3.4, the CS charges are 
calculated. Section 4 contains our numerical results. In 4.1 we have constructed the noncontractible 
loop (NCL) for the unstable solutions and have plotted the energy of the solution versus the angle 
parametrising this loop. In those cases where there is a nonvanishing CS charge, the energy is plotted 
also against the CS number. In 4.2 we have constructed the negative modes excited by the instability 
and have calculated the negative eigenvalues. Section 5 is devoted to a summary and discussion of 
our results. 

2 Should it turn out that the gauge group must be SO(n), n < D, it would only be possible to construct solutions 
subject to less than spherical symmetry. 
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2 The models and Chern-Simons densities 



In the first subsection 2.1 we define the static energy density functional resulting from the Lagrangian 
(1) in spacetime dimensions d — 6, 7, 8, and consider the topological lower bound on the energy for 
the d — 7 case. In the second subsection 2.2 we define the Chern-Simons densities for the models in 
spacetimes d = 6 and 8. 

Spacetime coordinates are labeled by Greek indices p,,u,..., and spacelike coordinates by Latin 
indices i, j,k, .... 



2.1 The models 

Since we are interested in static solutions, we will define the models on the d—1 dimensional Euclidean 
spacelike manifold. In the present work, we will restrict our considerations to the P = 2 model (1) 
used in [11]. Expressed in component form this is 

n * = -ik™ F » + £i T[F w (2) 

where 

Fijki = {Fiu,F k i]} 

is the 4-form curvature, and [jkl] implies cyclic symmetrisation. In the light of Derrick's scaling 
requirement, (2) is the simplest model that can support static finite energy solutions for spacetime 
dimensions up to d — 8. Beyond that terms with higher values p must be included in (1), but this is 
entirely unnecessary since all our conclusions from the present investigation hold qualitatively also 
for d > 9. 

The model is specified finally by the choice of the gauge group G, as well as its represent aion. 
We will mostly restrict ourselves to G — SO(d), but for odd d we will include the special case of 
G — SO(d - 1) as well. 

When d is even, i.e. d — 6 and 8 in our examples, there exists a Chern-Simons (CS) density in 
the 5 and 7 spacelike dimensions and the static solutions will display a sphaleron like unstability. As 
will be seen below, the nonvanishing of the CS density is what necessitates the choice of G — SO(d), 
with d — 6 and 8. These two examples will be the main focus of our attention. 

When d is odd, i.e. d — 7 in our examples, there exists a Chern-Pontryagin (CP) density in the 
6 spacelike dimensions. Provided that G and its representation is chosen suitably, this CP density is 
nonvanishing and the resulting static solution will be stable. As will be seen below, the existence of 
a stable soltion is what necessitates the choice of G — SO(6)(— SO(d — 1)) here. Otherwise, with 
G = SO (7) (= SO(d)), the static solution will not be stable. 

Before proceeding to state the CP and the CS densities for the even d systems, we give the 
topological lower bound on the energy density functional of the odd d, namely the d = 7 system. 
The inequality 

Tr {^JrlFij — ^^ L e i j k i mn F k i mr }j > (3) 

leads to 

> o -\/ r l r 2 £ ijklmn Tt FijFklF mn , (4) 
o 

stating the topological lower bound in terms of the 3rd CP desity. It follows that the spherically 
symmetric solution [16] to the model in d = 7 is topologically stable provided that the CP density 
on the right hand side does not vanish. This is the case for G = SO(d — 1) = 5*0(6), but only if the 
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gauge fields are in the chiral represenation of G. These static solutions in even spacelike dimensions, 
which obey instanton like boundary conditions, play the role of solitons. The topological inequality 
(3) cannot be saturated since the corresponding Bogomol'nyi equations are overdetermined. 



2.2 The Chern-Pontryagin and C hern- Simons densities 

The definition of Chern-Simons (CS) densities in d — 1 spacelike dimensions follows from that of the 
Chern-Pontryagin (CP) densities, in even spacetime dimensions, 



Qd — AT £ MiM2M3---Md^ r ^MiM2^2M3----^)td-iM<i ' ^ 6 and d 8, 
l*d 



(5) 



where N d is the appropriate normalisation factor for dimension d. Of course, for the purpose of 
definition (5), the signature is taken to be Euclidean. 

The CP density is a total divergence, which we denote generally as 



In the two cases of interest at hand, with d — 6 and d — 8, Qffl is given by 



(6) 
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From (7) and (8), the respective CS densities va-i are defined to be the /x = components of in 
5 and 7 Euclidean spacelike dimensions 
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(9) 



(10) 



3 Imposition of spherical symmetry 

Our choice of gauge group will be G = SO(d), except in d — 7 where we will consider and dispose 
of the special case G = 5*0(6) = SO(d — 1). Our choice of G = SO(d) for even d is similar to 
the choice for the Weinberg-Salam [17] (WS) and for the Bartnik-McKinnon [21, 22] sphalerons in 
d = 4. Furthermore, like in the latter case, we will employ the chiral representation G = SO±(d), 
noting that for the special case d = 4, G = SO±(A) = SU L / R (2). 

Since our main aim is to study the instability of the solutions to (2), we will employ the general 
spherically symmetric Ansatz parametrised by three radial functions u>i(r), u>2( r ) and w^(r). This 
follows from the axially symmetric Ansatz [20] in d dimensions where all the components of the gauge 
connection are taken to be independent of x and the component A = 0. Like the sphalerons [17, 22] 
in d — 4, our solutions are also parametrised by the function w\ while the functions (w 2 ,w 3 ) will be 
excited only in the directions of the instability. 
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For even d with G = SO(d), including d — 4, the imposition of spherical symmetry on the gauge 
field on (d — 1) dimensional Euclidean space results in 



Ai — Xj H -— (5ij — XiXj) Ej-J H -— XjXj E^J (11) 



where Il[f is (one of) the chiral representaion(s) of the algebra of SO(d), namely 

(d) ( i ± r^+i 



Eg' = (^f^Jly (12) 
if = —[Ti.r,-], (is) 



r^- being the spinor represention matrices of the algebra of SO(d) defined in terms of Tj, the gamma 
matrices in d dimensions. 

For odd d with G = SO(d), there are no chiral representations, so the Ansatz (11) holds only in a 
formal way, replacing the matrices by if. 

For odd d with G = SO(d), it is possible to employ the chiral representation of SO(d — 1). Here 
again, (11) holds formally, but now by replacing d with d — 1 in it. As noted above, in this case the 
solution (with w 2 = w 3 = everywhere) will describe a stable soliton. 

For the evaluation of the residual one dimensional energy density functional, only the algebraic 
properties of E^ and of enter the calculations, so the same result holds (up to an unimportant 
numerical factor 2) both for even and for odd d and with all G. For the evaluation of the CP and 
CS densities however, this distinction must be kept. 

3.1 Reduced one dimensional systems and gauge freedom 

Imposition of spherical symmetry on the (d — 1) dimensional SO(d) system results in the one dimen- 
sional energy density functional 



d 2-2! 4 



/ , w 2 w 3 \ 2 ( , W!W 3 \ 2 d-3 2 2 



d-2 

+ w 32 {d ~ 2){d - 3){d - 4) r "~ 8(1 - N2)2 



/ , . w 2 w 3 \ 2 ( , wxwzy d-5 2 J 



(14) 



in which we have used the shorthand notation 



I 2 1 2 

= W 1 + W 2 . 

For odd d = 7 and with G = 5*0(6), (14) holds as well, but the solutions to the field equations 
(to be presented in the following subsection) are quite different, with w\ = w 2 = everywhere. 

In the generic case, for the fields (11), the requirement of analyticity at the origin r = results 
in the asymptotic conditions 

wi(0) = l , w 2 (0) = , w 3 (0)=0, (15) 
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while in the asymptotic region r ^> 1 the requirement of finiteness of the energy results in the 
boundary condition 



lim \w\ 2 = lim (w 2 



W. 



2\2 



1. 



(16) 



(16) can be parametrised, like in [19], in terms of the familiar angle q as 

lim w\ = cos a , lim W2 = sing. (17) 

The conditions (16)-(17) can be understood better by displaying the gauge freedom of the energy 
density functional (14), under the action of the local U(l) transformation 

g[u{r)\ e U(l) , limcj(r) = , kw(r) = -q . (18) 
The action of this U(l) gauge transformation is given by 



w\ \ g[u>] 1 cos uj —sin a; \ I w\ 
W2 J \ sincj cos uo J \ W2 

w 3 1— > w 3 + ru; . 



(19) 
(20) 



3.2 The Euler— Lagrange equations 

The Euler-Lagrange equations of the SO(d) systems are, for variations of wi, W2 and w 3 in that 
order, 

ri {< + ^[{d-A)w' 1 + 2w' 2 w 3 + w 2 w' 3 } + ([(d- 3)(1 - |w| 2 ) - w^wx + (d - 5)w 2 w 3 ) } 
+ ||(d-3)(d-4)(l-M 2 ). 

{(1 - H 2 ) (w? + J[(d - 8)w[ + 2w' 2 w 3 + w 2 w' 3 } + ^ ([(d - 5)(1 - H 2 ) - w 2 ]wi + (d - 9)w 2 w 3 )^J 



+ 2 



w[ + 



+ w 



., n -pV ( 1 . ^ 2 W 3 ^ 
lUl - 2 ( |w| ) \w 1 H J 



= 



(21) 



+ 



{wZ + ^[(d-4)w' 2 -2w' 1 w 3 -w 1 w' a \ + ^ ([(d-3)(l - \w\ 2 ) - wl]w 2 - (d - 5) Wl w 3 ) J 



3r 2 
4^ 



(d-3)(d-4) (1 - \w\ 



|(1 - |w| 2 ) + h(d- 8)w' 2 - 2^3 - + ([(d - 5)(1 - \w\ 2 ) - wl]w 2 - (d - 9)w!W 3 ) 



+ 2 



/ , ^2^3 

w 1 H 



3r 2 



) +{w 2 -—) w 2 -2(\w\ ) [w 2 -— — J I =0 



(22) 



r 1 + — (d-3)(d-4)(l-H 2 ) 2 



(wiw' 2 — w 2 w[) w 3 \w 



■n\2 







(23) 



The last equation, (23), is satisfied by the solutions of (21) and (22), as can be seen straighforwardly 
by differentiating (23) and identifying it with the difference of w 2 times (21) and w\ times (22). Thus, 
exploiting the gauge freedom (19)-(20), we choose the radial gauge with w 3 = 0. We then proceed 
to solve the two equations (21) and (22) for the two functions W\ and w 2 . 

We note here that in the w 3 = gauge, equations (21) and (22) are symmetric in the functions 
w 1 and w 2 as in the case of the Bartnik-McKinnon sphaleron [21, 22], and in contrast to the WS 
case [19] in which this symmetry is absent due to the presence of the complex doublet Higgs field. 
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3.2.1 SO±(d— 1) solitons in d dimensions 



Before proceeding to consider the sphaleron solutions, we dispose of the stable soliton solutions of 
the models in odd d spacetime dimensions when the gauge connection takes its values in the chiral 
representation of SO(d—l). We restrict this demonstration to the model in d = 7, with G = SO±(6). 

As a practical illustration of the topological lower bound (4), using the notation K-y = \J (d — 2)ti/2, 
«2 = 3a/ (d — 2){d — 3)(d — 4)r 2 , consider the two inequalities 



> (24) 



> . (25) 
Adding (24) and (25) we have 



H d [ Wl , w 2 , w 3 ] = H d [w, 0, 0] > 2k,k 2 \J ^ + y ^ j r d - 7 (1 - w 2 ) w' , (26) 

where Hd[w, 0, 0] is given by the energy density functional (14) with w± — w and w 2 = w 3 = 0. It is 
now obvious that the right hand side of the inequality (26) is a total derivative only for the spacetime 
dimension d — 7 at hand, and that it is the residual CP density after the imposition of spherical 
symmetry on the G = SO±(d — 1) = SO±(Q) model, with G is in the chiral representation. 



d •"» . 2 \ 2 



1 2 



Kl r 2 w ' - K2r 2 J — (i _ w ) 



d-3 



K 2 r 2 (1 — w )w' — K,±r 2 y - - (1 — w ) 



3.3 Reduced C hern- Simons densities and non— contract ible loop 

The reduced one dimensional Chern-Simons (CS) densities of static gauge fields (11) in (even) d 
spacetime dimensions, with gauge groups G = SO(d) in the chiral representation SO±(d), take on 
nonvanishing values 3 . The particular examples considered concretely here are those in d = 6 and 
d — 8, given by (9) and (10). Subjecting these to spherical symmetry according to the Ansatz (11), 
we find 

5 

n 5 



n 7 = 



- A { [2(1 - Wl ) + 3(1 - H 2 )] [(iui - IX - w 2 w[] + *w 3 (1 - H 2 ) 2 } (27) 
. _A^{ [5((1 - ^) 2 + W 2 ) 2 + 6(1 - «; 1 )(3(1 - w x ) + 4(1 - \w\ 2 )}] [( Wl - l)w' 2 - w 2 w[] 

+ ^(l_ H ')3}. (28) 



In what follows a most important role will be played by the non-contractible loop (NCL) of 
configurations displaying the instability of the sphalerons. In contrast to the case of the WS model 
[19], and similar to the Barnik-McKinnon sphaleron [21, 22]. This is a result of the (wi, w 2 ) symmetry 
in the field equations (21)-(22). In the former case where the doublet Higgs field removes this 
(wi,w 2 ) symmetry, conditional solutions with fixed CS number Qcs can be constructed by solving 
the equations of motion that extremise the energy density functional plus a Lagrange multiplier times 

3 Note that if the representations of G adopted were not the chiral ones given by the spin matrices (12), but rather 
those given by (13), then the resulting CS densities would vanish. 



8 



the CS density. In the that case [19], the minimal finite energy path versus Qcs can be constructed 
concretely. The situation here is similar, rather, to the second case of the EYM sphaleron [21, 22], 
where the analysis of instability is carried out exclusively by employing a NCL. 
The actual sphaleron solutions are parametrised by the functions 

W! = f(r) , w 2 = , w 3 = (29) 

which will be constructed numerically in the next section, consistently with the asymptotic conditions 
(15) and (16), the function /(r) in (29) satisfies the asymptotics 

lim/(r) = l , Hm/(r) = -l. (30) 

Following [23] we adopt the NCL configurations 

1 1 

wi = -(l + cosg) + -(l-cosg)/ (31) 

w 2 = l(sing)(l-/) (32) 
w 3 = 0, (33) 

which by virtue of (30) ensures that both (15) and (17) are satisfied all along the NCL. 
3.4 Calculation of the C hern— Simons number 

In the case of the WS model [17, 18, 19, 22, 21], the CS number Qcs has the physical interpretation 
of Baryon number. Here, it merely is a convenient topological charge characterising the sphaleron. 

The topological CP charges are are given by the d dimensional volume integrals of d^Q^ . Here we 
are interested in the k = 3 and k = 4 examples for which are given by (7) and by (8) respectively, 
but now we consider static fields in a d dimensional Minkowskian spacetime. The topologial charge 
in this context is referred to [18] as the CS charge Qcs- 

Adapting the arguments of [18] to our cases, notably assuming that Qcs — at t — — oo, the 
topological charges are given by 

*» = m (/'I * jL ^ dS ' + L dd " x • (34) 

The Chern-Simons (CS) densities = Vd-\ are displayed by (9)-(10). The task at hand is to 
evaluate the integrals (34) for our two sphaleron solutions. 

In the w 3 = gauge in which the sphaleron solutions are constructed, the surface integral term 
in (34) does not vanish since the gauge potential decays with power r _1 as r — > oo, as seen from (11) 
and (29)- (30). It is convenient to evaluate (34) in a gauge in which the surface integral vanishes, 
such that 

= m L dd ' lx ^ = w) L Vi - 1 dd ' lx • (35) 

which for the spherically symmetric fields (11) reduces to the one dimensional integral 

r-r=oo 

Qcs = / n d _i dr (36) 

Jr=0 

with n d _i(r) given by (27)- (28). 
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But substituting the NCL configuration (31)-(33) in the expressions for CS densities (27) and 
(28), results in the vanishing of and of n-j. This means that the nonvanishing contribution to the 
integral (34) must come from the surface integral term, which is not convenient to evaluate since the 
solutions and the NCL configurations at our disposal are time independent. But by definition (34) 
is gauge invariant, so it should be evaluated in a gauge in which the surface integral in (34) does not 
contribute. To this end, following [18], we subject (31)-(33) to the gauge transformation (19)-(20), 
such that 

lim u)(r) = , lim u(r) = —q . (37) 

With these boundary conditions on the U(l) gauge group parameter, and (30) for the sphaleron 
profile function f(r), we find 

lim u w\ = 1 , lim W W2 = , lim w -uJ 3 = , (38) 

r^oo r^oo r^oo 

which results in the connection (11) decaying faster than r~ x at infinity. As desired the surface 
integral term in (34) now vanishes. The density n d -i in (36) does not vanish, and can be evaluated 
to give the CS number. 

Subjecting the NCL configuration [w\, W2, ws) given by (31)-(33) to the gauge transforamtion 
(19)-(20) with (37), and substituting the resulting set of functions ("wji," w 2 , w w^) into the densities 
(27) and (28), we obtain total derivatives in the variable r. After integration, we have the respective 
CS numbers 4 

Qcs = ^ (3g-4sing + isin2g) (39) 
21 

Q cs = — (30g-45sing + 9sin2g-sin3g) (40) 

207T 



4 Numerical results 

4.1 The sphaleron on the NCL 

We have solved the system (21)-(22) numerically for d = 6,7, 8, i.e. for the gauge group G = SO(d). 
(This excludes the case of odd d with G = SO(d — 1) with topological lower bound (26)). The 
coupling constants Ti,r 2 can be chosen arbitrarily by choosing an appropriate scale of the mass M, 
defined as the integral of H d 

M = J H d dr , 

and of the radial variable r in (14). We have chosen t\ — r 2 = 1. However, for practical reasons, 
the numerical values for the masses and for the negative modes given below will be given in units 
2( d ~ 2 ^ 2 S d - 2 , where S d ^ 2 denotes the surface of the sphere in the d-dimensional space-time. 

The three profiles for the function Wi(r) are presented in Figure 1, and the figure reveals that the 
dependence of the profile on d is rather weak. 

In the units choosen, the masses of these solutions are respectively M(d — 6) ~ 4.67, M(d — 7) ~ 
11.3, M(d = 8) ~ 36.8. From the numerical profile for w±, the different configurations of the path 
(31)-(33) can be constructed; their energies can further be computed as functions of the parameter 
q. The energy is plotted as a function of the parameter q on Figure 2 for d = 6, 7, 8. As expected, 
the figure shows that the configurations on the path have finite energy and that the energy increases 

4 In d = 4, n 3 = ± [(wi - l)w' 2 - w 2 w[ + *f(l- \w\ 2 )] , yielding the familiar [19] result Q cs = SL if L£ . 
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Figure 1: The profiles of the function w\ for the sphaleron solutions in d — 6, 7, 8, which are practically 
unchanged when the function w<i is perturbed by 5w2- The profiles of the negative modes 5w2 are 
displayed and the corresponding negative eigenvalues u 2 listed. 
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Figure 2: The energy along the path extrapolating between the vacuum and the sphaleron is plotted 
as a function of 6 = q/n for d=6,7,8. 
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monotonically from the vacuum (q = 0) to the classical solution (q — 7r), demonstrating that the 
classical solution is indeed a sphaleron. 

For the cases with even d and nonvanishing CS densities (27)- (28), the energy is plotted also 
against the CS charges (39)- (40) respectively, on Figure 3. 

4.2 Negative modes 

In the second subsection we have carried out an infinitesimal stability analysis by constructing the 
negative modes and evaluating their negative eigenvalues. We also track, within the limits of validity 
of this infinitesimal analysis, the growth of the CS charge as a function of the w\ and the function w 2 
excited along the direction of instability. For this purpose we consider a perturbation of the classical 
solution constructed in the previous subsection : 

w l = wi + e iu3t 5wi , w 2 = w 2 + e iuJt Sw 2 , w 3 = , (41) 

choosing again to work in the gauge w 3 = 0. Inserting the perturbed solution into the equations 
of motion and retaining only the linear terms in the fluctuation, we got a system of two decoupled 
Sturm-Liouville equations in 5w\ and Sw 2 . The construction of the normal modes is equivalent to 
finding the normalisable solutions of these equations, which is a problem beyond the scope of this 
paper. We limit our analysis to the construction of the mode of lowest eigenvalue. The results of 
[24] strongly suggest that the main mode of instability of our solution should appear in the sector 
5wi = 0. This turns out to be the case. Technically the negative mode can be constructed by 
minimising the quadratic form 

5E = 2^ 2 S d J-^- J drr d ~ 4 (h{5w' 2 f - I 2 {5w 2 f) (42) 

with 

h = fa + 3r 2 (d - 3)(d - 4)1(1 - w\f) (43) 

h = {d- 3) (1 ~ Wl2) (n + 3r 2 (d - 4)l((d - 5)(1 - wlf + 2r 2 K) 2 )) . (44) 

In (43) -(44), the function w\ represents the classical sphaleron solution. Inspecting the form of this 
variational problem we see that it leads to a Sturm-Liouville equation with a potential given by the 
function —I 2 . This corresponds to a potential well and allowing the existence of a negative mode. 
We were able to construct numerically one normalised negative mode Sw 2 in each case d = 6,7, 8, 
whose profiles are displayed in Figure 1, together with the profiles of the corresponding classical 
solution w±. The fact that 5w 2 presents no node strongly suggests that our solution corresponds 
to the eigenmode of lowest eigenvalue. These eigenvalues, indeed, appear to be negative and were 
evaluated to be — u 2 = 3.88, 17.4, 84.6 respectively for d=6,7,8, again using the same scale as before. 

Finally we evaluate the CS charge Qcs for d = 6 case, for the configurations of the form w\ = w\ 
, w 2 = e5w 2 , e being an infinitesimal parameter. We find 

Qcs(wi, w 2 ) = Qcs{q = tt) - 1.54e (45) 

where Qcsil = vr) is the CS charge of the sphaleron, which for d — 6 is evaluated from (39), as 
Qcs{q = vr) = 22.5. This shows that the CS charge varies linearly with the parameter e when the 
sphaleron is perturbed in the direction of its unstable mode. 
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Figure 3: The energy along the path extrapolating between the vacuum and the sphaleron is plotted 
as a function of q for d=6,8. 
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5 Summary and Conclusions 



We have studied static finite energy solutions to Yang-Mills systems in higher dimensions. The 
interest in these solutions is that they are gravity decoupling limits of the fully gravitating Einstein- 
Yang-Mills solutions in higher dimensions. In turn, gravitating YM systems in higher dimensnions 
are important field theoretic models arising in the study of Dp-branes in spacetime dimensions 
larger than 4, in the context of the low energy effective action of superstring theory and gauged 
supergravities. Some of these solutions have sphaleron type instabilities, which must be studied. In 
this respect the higher dimensional gravitating YM systems are similar to the gravitating YM model 
in 4 spacetime dimensions, whose sphaleron like instabilities were studied in detail in [21, 22]. What 
is different in the higher dimensional cases at hand, versus the 4 spacetime dimensional case [21, 22], 
is that unlike for the latter, here we have gravity decoupling limits, and, the sphaleron nature of the 
gravitating solutions is essentially identical to their flat space counterparts. The sphaleron analysis 
being appreciably simpler in the flat case, we have chosen to work with those models. 

We have restricted our studies to that of spherically symmetric solutions only. This has neces- 
sitated certain, rather limited, choices of the gauge groups G. These gauge groups turn out to be 
SO(d— 1) and SO(d) in odd spacetime dimensons d, and SO(d) only in even d. The gauge connections 
take their values in the spinor reprentations in terms of Dirac matrices, and whenever G = SO(N) is 
defined for even N, we have employed the chiral representations of SO(N). All solutions considered 
are evaluated numerically and they satisfy the second order Euler-Lagrange equations rather than 
first order self duality equations, since the YM systems in question are not scale invariant. 

We have shown that in odd spacetime dimensions d, the solutions to SO(d — 1) YM systems are 
topologically stable if the representaion of G employed is chiral, i.e. that they are solitons which 
are stabilised by the Chern-Pontryagin (CP) topological charge. These are the only models which 
support stable solitons. Should the representations of the SO(d— 1) algebra employed not be taken to 
be the chiral ones, the solutions will be unstable. They are unstable also for the choice of G = SO(d). 

In even spacetime dimensions, the solutions turn out to be always unstable, due essentially to the 
absence of a CP topological charge in the odd spacelike dimensions. For G = SO(d), by employing 
the chiral representations, we have calculated the Chern-Simons (CS) charges, and have plotted 
the energy of the noncontractible loops versus the CS charge Qcs, highlighting the nature of these 
solutions as types of sphalerons. 

For all the unstable solutions, namely to the SO(d) models in all d dimensional spacetimes, we 
have constructed the negative modes of the corresponding fluctuation equations and calculated their 
negative eigenvalues. For the SO±(6) solutions in d = 6 in particular, we have also calculated Qcs 
perturbatively, showing that in the region of validity the Qcs increases linearly in the pertubation 
parameter. 

We have restricted our study to the simplest model (1), with p = 1 and p = 2 terms only, 
in spacetime dimensions d = 6, 7, 8. This limited choice is sufficient to illustrate the questions of 
stability and instability in the generic cases. The properties demonstrated in the p = 2 examples 
studied repeat themselves in every 4p dimensions. Thus for example for p = 3 models, in spacetime 
dimensions d — 6, 8, 10, 12 the solutions will be unstable sphalerons, while those in d — 7, 9, 11 will be 
stable, being stabilised by the 3rd, 4th and 5th Pontryagin charges, respectively. If the d = 7 model 
chosen in that case does not contain the p = 2 (in addition to the obligatory p = 3 term dictated 
by Derrick)term in the YM hierarchy, then there will be no gravity decoupling limit. In that case 
the sphaleron analysis must be carried out in the fully gravitating model, following the lines of the 
analysis in [21, 22] for d — 4. The same holds true in the case of the p = 2 model studied here, d — 5. 
In that case too there exists no flat space limit so that the sphaleron analysis must again be carried 
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out as in p — 1, d — 4 [22, 21]. The latter analysis falls outside the scope of the present work and is 
deferred to a future study. 
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